Journal of Research of the National Bureau of Standards 



Vol 48, No. 6, June 1952 



Research Paper 2332 



Two Applications of Group Characters to the Solution 
of Boundary- Value Problems' 



E. StiefeP 




It is shown that the numerical work involved in solving a boundary-vame or eigen- 
value problem by finite difference methods in a domain with many symmetries can often 
be reduced by applying the theory of group characters to the group of symmetries of the 
domain . 

In the second part the author considers the problem of solving Am — in a cube when 
the prescribed boundary values are invariant under the group of rotations of the cube. 
In the series representation of the solution in terms of harmonic polynomials only a subset 
of these polynomials actually occurs, and the theory of group characters facilitates con- 
siderably the determination of this subset. 

It is very well known that representation of groups 
may be a useful tool in discussing various eigenvalue 
problems in quantum mechanics. This is due to the 
fact that boundary-value problems of Schroedmger's 
equation are to be solved mostly in the whole space 
or inside a sphere as fundamental region and therefore 
the high symmetry of those domains can be used. 

It is the purpose of this paper to show that group 
characters sometimes are appropriate also to sim- 
plify the numerical computation of boundary-value 
or eigenvalue problems in more general-shaped do- 
mains, in particular if the problem has been trans- 
lated into the language of difference calculus. It 
will be sufficient to explain this discussing simple 
examples. 

Let, for instance, A be any plane region and D a 
differential operator invariant under Euclidean mo- 
tions or reflections of the plane. More exactly, if 
J{x,y) is a function in our plane and M a Euclidean 
rotation, translation, or reflection, then 



D(Mf) = M{Df), 



(1) 



which is to say that D must commute with any 
motion M. The Laplace operator A and the operator 
AA of elasticity problems have this property. 

In order to solve any given boundary-value prob- 
lem with respect to the region A and the operator 
D, we may use the difference technique, introducing 
a square lattice L in A and replacing / by a lattice 
function (defined only in the lattice points) and D 
by a difference operator, which approximates the 
given differential operator. In the case of a A-prob- 
lem, the corresponding dift'erence operator may be 
given by figure 1. Relation (1) will be true again in 
the difference case, provided that by Mis understood 
any motion or reflection transforming the lattice 
into itself. 

Suppose now^ that the given region A is invariant 
under a group G of transformations of this type. If 
A has no symmetry at all, we may use Schwarz's 
alternating process dividing A into overlapping sub- 
regions, each of them having symmetry properties. 

1 This work was performed on a National Bureau of Standards contract with 
the University of California, Los Angeles, Calif., and was sponsored (in part) 
by the Office of Naval Research. 

2 University of California, Los Angeles, Calif., and Eidgenossische Tech- 
nische Hochschule, Zurich, Switzerland. 
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Figure 1. 

It is our goal to show how the group G can be used in 
order to simplify the solution of the given boundary- 
value problem in the region A. 

Take, for instance, the square region A of figure 2. 
The group G is in this case the symmetry group of 
the square of order 8 having the following similarity 
classes : 



Class 


Elements 


Number 

of 
elements 


[1] 

[2] 

[3] 

[4] 

[5] 


Identity 

Rotations ±90°__ __ __ _ 


1 
2 

1 

2 
2 


Rotation 180° 


Reflection in horizontal or 
vertical line 


Reflections in diagonal lines _ 
Total 


8 



(2) 



The given boundary-value problem is equivalent to 
a system of nine linear equations for the values of the 
wanted function m the nine interior lattice points. 
Any lattice function / may be visualized as a vector 
in a nine- dimensional vector space S. Therefore, 
applying the operations of the table to the lattice 
functions, we get a linear representation R of the 
basic group G in S. 
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Figure 2. 

It is not too hard to show that its character x is 
given by the following rule: The value x(M) of x for 
an element Al of G is equal to the number of lattice 
points fixed under the operation .1/. (Use as coordi- 
nate system in S the lattice^ functions not equal to 
only in one sing'le lattice point; they are })erniuted l)y 
M.) We want to decompose^ the represcMitation Ii 
into irrculucible conipoiUMits: 



H = C,I^,+C,fi2 + 



+cJh, 



(3) 



Where /?i, /i^, . . ., R^ are the irreducible^ re])re- 
sentations of G. (There are exactly five because G 
has five similarity classes). According to the oi'lho- 
gonality theorems of representation theory this has 
to be done working on the table of characters of G 
ns follows : 



Class 

No. ofelements_ 


[11 [2] [31 [41 [5] 
12 12 2 


Xi 

X2 

Xa 
X4 
X5 


11111 

1 1 1-1-1 
1-1 11-1 
1-1 1-1 1 

2 0-200 


X 


V) 113 3 



(4) 



In the first line we have listed the similarity classes 
with their numbers of elements and in the following 
lines the values of the characters of the irreducible 
representations in those classes. They are recorded 
in many textbooks of representation theory.^ The 
last line is the (diaracter x of our representation R, 
computed according to the rule previously mentioned. 
The first number in each line is the degree of the 
representation (dimension of the vector space trans- 
formed by the representation). According to a 
general formula of character theoiy the coefficient 
Ci in (^) is 



Ci=ij:x{M)xi{M), 



(5) 



with the group element M running through the whole 
group G. By this we have for instance from our (able 

C4=--i(9-l -2-1 + M -2.3 + 2.3) = !. 

The multiplicities c, are listed in the column to the 
right of the table. Finally, we find 

H=U\ + R^ + h\ + 2lh. (6) 

Therefoi-e, the vector space S splits into subspaces 



Si, Si, 



S:i', S^; Ss, S^, 



a) 



each of them being invariant under R and being 
transformed by the ir'reducihle representation given 
by its subscript. 

Take now a set of subspaces having the same sub- 
script (i. e., being transformed by the same repre- 
sentation) and introduce base vectors in each space 
of the set. For instance, S5 and S^^ being both 
two-dimensional, may be spanned by the base vectors 



S,: 



('52 



(8) 



3 A. Speiser, Die Tlieorie der (Tiiippen von endlicher Ordnung, 3d ed., p. 179- 
180 (Springer, Berlin; and Dover Pub). 



Generally speaking, each set of (Hjually transformed 
subspaces yields a rectangular arrangement of base 
vectors by listing the base vectors of each subspace 
in a line. The length of the rectangle is the degree 
of the irreducible representation under consideration 
and its height the multiplicity of this representation 
as a component of R. It is, of course, possible to 
choose the base vectors in such a way that each line 
of the rectangle is ti'ansfornied exactly in the same 
way. 

It follows then from a well-known generalization ^ 
of Schur's lemma that each linear operator in S com- 
mutable with the transformations of I\ transforms 
the vectors of each column in our rectangle among 
themselves and in each column in the same way. 
Taking into account the relation (1), we see that the 
dift'erence operator of the given boundary-value 
problem has this property. From this follows imme- 
diately that the original system of nine liiu^ar equa- 
tions with nine unknowns splits into partial systems. 
The rectangle (8), for instance, yields two systems of 
two equations each, having the same coefficient 
matrices. From the first and the c column of the 
table (4) we get therefore the following final result. 
The system of nine equations splits into one system of 
three simultaneous equations, one single equation, 
one single equation, and two systems of two simul- 
taneous equations. 

If we have to deal with an eigenvalue problem in- 
stead of a l)oun(larv-value problem, the analogous 
statement is, that the wanted characteristic poly- 
nomial is the product of a factor of third degree, two 
linear factors and the square of a quadratic factor. 
In order to carry out the decomposition of the system 
of nine equations explicitly, it is, of course, necessaiy 

< B. L. van der Waerden. Die gruppentheoretische Methode in der Quanten- 
mecluinik, sect. i:i, p. 47-50 (Springer, Berlin, 1932). 
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to establish the lattice functions, which build the 
base vectors of the different rectangles (8). The 
three basic lattice functions corresponding to the 
first irreducible representation /?i are simply three 
linear independent functions invariant under all 
motions of the square. As another example we give 
the values of the lattice functions (8) in the nine 
lattice points. 





^51= —1 1 , ^52 = 



-1 1 

eL= 0, el,-- 

-1 1 



Taking as /} the operator of figure 1 , we get 
De',,= 2^51-4^;,. 






1 

















— 1 





1 





1 











1 





-1 



(9) 



(10) 



The functions e^o, ^52 are transformed in the same way. 
The eigenvalue problem of the operator D with 
vanishing boundary values gives the equation 



Df+\f = 0. 



(11 



Let us put /=a€5i + 6651. Then (10) yields the eigen- 
value problem 

(X-4)a + 2ft-0 

(12) 
a+{\-4:)b = 0. 



And from this follow the eigenvalues 
X = 4±V2 



(13) 



of the operator D, each of them having at least the 
multiplicity 2. 

It is not too hard to establish for a lattice in a 
given region A the rectangles of base vectors when 
we know the characters and the irreducible repre- 
sentations of the symmetry group oi A. It is con- 
venient to choose all those vectors (their total number 
is equal to the number of lattice points) orthogonal 
to each other. This is possible because the re- 
presentation R is orthogonal, and therefore the ir- 
reducible subspaces of a§ are orthogonal. 

At the Institute of Applied Mathematics in Zurich 
(Switzerland) we used the methods of this paper in 
order to solve the boundary-value problem of 
Airy's elasticity function inside a dam. As subre- 
gions in Schwarz's alternating process, we chose 
squares of 16 lattice points. 

In his book about eigenvalue problems, L. Collatz ^ 
introduces the notion of '^eigenvalues of a graph" 
and discusses an example, where the graph is built 
by the segments joining two midpoints of 
edges of a cube. This problem ma}" also be treated 

5 Jj. Collatz, Eigenwortaufgaben mit technischen Anwondungen, p. ;:82 
(Leipzig. 1949). 



by characters, using the group of all rotations of 
the cube. 

In this section another application of character 
theory is discussed, related to the methods of S. 
Bergman and M. Picone for the solution of 
boundary-value problems. For our purposes those 
methods may be characterized by the following 
procedure. At first a complete set of particular 
solutions of the given partial differential equation 
is constructed and then the solution of the boundary- 
value problem is expanded into a series having as 
terms those particular integrals. Let us take, for 
instance, Laplace's equation Aii = m three-dimen- 
sional space. The particular solutions may be in 
this case the harmonic polynomials of degree n 



Y^-\x,y,z)-- 



,y>«.g^^t^.p^-)(cos (9), 



n=0,l,2, 



(14) 



m running from ( — n) to {-\-n). On the right side, 
r, 6, </) are polar coordinates, P^'"^ is a generalized 
Legendre polynomial. Suppose now^ that the fun- 
damental region A of the boundary-value problem 
is a cube having its center at the origin and suppose 
furthermore that the given boundary values are 
invariant under the group G of the 24 rotations of 
the cube. It is obvious that under those circum- 
stances only harmonic polynomials appear in the 
Bergman-Picone expansion, w^hich have the same 
symmetry properties. In other words, we have 
the problem to establish all harmonic polynomials 
invariant under the group (7 of rotations of the cube.^ 
The theory of characters gives general methods 
in order to solve problems of this type. G has the 
following similarity classes : 



Class 


Typical element 


Number 
of ele- 
ments 


[1]-- 
[2... 
[3__. 

[5]-.- 


Identity 

Rotation ±90° 


1 
6 
3 

6 

8 


Rotation 180° (z-Sixis) 

Rotation 180° (axis: z = 0,x=y)_. 
Rotation ± 120° (axis = diagonal) _ 

Total 


24 





(15) 



By any rotation of the cube the (2n+l) harmonic 
polynomials of a given degree n are transformed 
linearly among themselves and yield therefore a 
representation R,, of degree (2/1 + 1) of G. Let M 
be any rotation of G and Xn{M) the character of Rn- 
The number of linear independent harmonic poly- 
nomials of degree n invariant under G is equal to 
the multiplicity c^ of the unit representation in the 
decomposition of Rn into irreducible components. 
(The unit representation of a group maps every 

6 This question has been posed and answered by G. Polya and B. Mever, 
Sur les svmetries des fonctions spheriques de Laplace, Compt. rend. 328, 28-30 
(1949); Snr les fonctions spheriaues ae Lay lace de symetrie cristallographique 
donnee, Compt. rend. 238, (1949\ 
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group eleineiit into the one-row matrix^l). This 
multiplicity is given by the character formula 



M running through the group. It remains to com- 
pute the character Xn{M). In order to do this 
we may restrict ourselves to rotations M around 
the s-axis. Let /z be the rotation angle. Taking 
into account (14), it follows that the (2n+l) har- 
monic polynomials of degree n are transformed by 
M in the following way : 

Y ir' -^e'^'^Y !t\ (17) 

The matrix of the transfomiation is pure diagonal 
with the trace 



+n 



Xn{M)=j:e'^^''= 



sin(2n+l)| 



(18) 



Trom (16) and the table (15) follows now the result: 



Cn=^{{2n 



+ l) + 6 



sin (2??.+ 1)45° 



+ 9 sin(2ri.+ l)90° + 8 



sin 45° 
sin (2n+ 1)60° 



sin 60° 



(19) 



« 



The final result may be expressed in i\\r following 
terms. Write n in the form of a mul(i|)l(' of 12 and 
tlie corresponding remainder 



Then 



n=l2'k+r, 0</<12. 

Cn=k+t{r), 
where t(r) is given by the table 



(20) 
(21) 



r 





1 2 3 


4 5 6 7 8 


9 


10 


11 


t{r) 


1 





10 10 1 


1 


1 






(22) 



The polynomials themselves can be built by choosing 
any harmonic polynomial of degrcn^ n, applying the 
24 rotations of the cube and adding up the 24 result- 
ing polynomials. 

A more general problem would be to establish the 
harmonic polynomials Ixung transformed l)y G 
according to a given irre(lucil)le I'epresentation of G. 
This can be used to simplify the solution of a bound- 
ary-value problem in which the given bounchiry- 
values have not the symmetry of the cube. The 
basic idea is, of course, the same as outlined previously 
in this paper. 



Los Angeles, January 29, 1952. 
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